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We present imaging observations at 1.3 milhmeters of the debris disk sur- 
rounding the nearby M-type flare star AU Mic with beam size 3" (30 AU) from 
the Submilhmeter Array. These data reveal a belt of thermal dust emission 
surrounding the star with the same edge-on geometry as the more extended scat- 
tered light disk detected at optical wavelengths. Simple modeling indicates a 
central radius of ~ 35 AU for the emission belt. This location is consistent 
with the reservoir of planetesimals previously invoked to explain the shape of the 
scattered light surface brightness profile through size- dependent dust dynamics. 
The identification of this belt further strengthens the kinship between the debris 
disks around AU Mic and its more massive sister star f3 Pic, members of the 
same ~ 10 Myr-old moving group. 

Subject headings: circumstellar matter — planet-disk interactions — stars: in- 
dividual (AU Microscopii) — submillimeter: planetary systems 
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1. Introduction 



van Leeuwen 



2007 ) MlVe flare star flTorres et al. 



^AU Mic is a nearby (9.91 ± 0.10 pc 

20061 ) with the young age of 12^*^4 Myr ( Zuckerman et al.ll200ll ). a key epoch in the forma- 



tion of planetary systems. The star is surrounded by a nearly edge-on circumstellar disk 
discovere d in coronographi c images of scattered starlight that extends to a radius of at least 
210 AU (jKalas et al.ll2004l ). Since the small grains responsible for scattering should be re- 
moved rapidly by stellar pressure forces, the disk is thought to consist of the collisional debris 
of unseen planetesimals experiencing ongoing impacts (for recent reviews of the debris disk 
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phenomenon, see IWyattI l2008l : lKrivovll2010l ). Because of its proximity, the AU Mic debris 
disk has become one of the best studied examples at optical and infrared wavelengths, in- 
cluding detailed iinagery of both inten sity and polarization from the Hubble Space Telescope 
ferist et allboosi ICraham et al.lboOTl ). 



Many observational properties of the AU Mic debris disk bear striking similarities to the 
archetype debris disk surrounding /3 Pic, also viewed nearly edge-on and located in the same 
young moving group. In particular, the midplane optical surface brightness profiles of these 
disks are remarkably alike, with shallow inner slopes (~ l/r^~^) that steepen substantially 
(^ l/r^~^) in the outer regions, near 35 AU and 100 AU fo r AU Mic and Pic, resp ectively 



flKrist et al.ll2005l : Fitzgerald et all 120071 : iHeap et al.ll2000l : iGohmowski et all 120061 ) . These 
similarities inspired the development of a unified framework for debris disks based on a 
l ocalized belt of planetes i mals, or "birth ring" , th at produces dust in a collisional cascade 
(IStrubbe fc Chiangll2006l : lAugereau fc BeustI 120061 ) . The smallest grains are blown out from 
the belt by radiation and winds from the star, larger grains are launched into eccentric orbits 
with the same periastron as the belt, and the largest grains, which are minimally affected 
by non-gravitational forces, do not travel far from the belt before being ground down. The 
grains are therefore segregated according to their size, which gives rise to the characteristic 
scatter ed light profile. For these dy namics to prevail, the disks must be substantially free 
of gas (IThebault &: Augereaul 120051) . a property con firmed by measurements at ultraviolet 



wavelengths (JRoberge et al. 



2005 



France et al.l 120071 ). 



A "birth ring" , if present, would remain hidden in optical and inf rared images dominated 
by small grains that populate an extended halo (e.g.. lSu et al.ll2005l ). By contrast, observa- 
tions at (sub)millimeter wavelengths highlight thermal emission from the largest grains and 
hence the location of the dust-producing planetesimals (I WyattI 120061). To date , the b est case 
for a "birth ring" comes from millimeter observations of (3 Pic (jWilner et al.ll201ll ). While 
the optical disk of (3 Pic extends more than 1000 AU from the star, the millimeter imaging 
reveals a much more compact belt of emission at ~ 95 AU radius. This millimeter emission 
belt coincides with the change in the slope of the optical surface brightness, which in the 
models marks the outer extent of the colliding planetesimals. 

Previous (sub) millimeter- wave observations of AU Mic did not have sufficient angular 
resolution to reveal much about the emission morphology. A detection at 850 /im using 
JCMT/SCUBA (14" beam) in dicated a reserv oir of cold dust with mass ~ 0.01 M®, but 
did not resolve any structure (JLiu et al.ll2004J ). Subsequent observations at 350 /zm using 
CSO/SHARC II (10" beam) margin ally resolved an or ientation compatible with the scattered 
light, but were otherwise limited ( IChen et al.l |2005[ ). In this Letter, we present imaging 
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observations of AU Mic at 1.3 millimeters using the Submillimeter Array (SMAJ^that obtain 
3" resolution and provide evidence for a planetesimal belt. 



Observations 



We observed AU Mic with the SMA (IHo et all 120041 ) on Mauna Kea, Hawaii at 1.3 mil- 



limeters wavelength using the compact and extended configurations of the array. Table [T] 
provides basic information about these observations, including the observing dates, base- 
line lengths, and atmospheric opacities. The weather conditions were good for all of these 
observations and best for the two compact configuration tracks when only 6 of the 8 array 
antennas were available. The phase center was a = 20''45"'09!53, S = -31°20'27?2 (J2000), 
offset by ~ 5'.'3 from the location of the star. The ~ 54" (FWHM) field of view is set by the 
primary beam size of the 6 meter diameter array antennas. The total bandwidth available 
was 8 GHz derived from two sidebands spanning ±4 to 8 GHz from the LO frequency. Time 
dependent complex gains were calibrated using observations of two quasars, J2101-295 (3.9 
degrees away) and J1924-292 (17.4 degrees away), interleaved with observations of AU Mic 
in a 15 minute cycle. The passband shape for each track was calibrated using available bright 
sources, mainly J1924-292. The absolute flux scale was set with an accuracy of ~ 10% using 
observations of Callisto or Ganymede in each track. All of the calibration steps were per- 
formed using the IDL based MIR software, and imaging and deconvolution were done with 
standard routines in the MIRIAD package. We made a series of images with a wide range of 
visibility weighting schemes to explore compromises between higher angular resolution and 
better surface brightness sensitivity. 



3. Results and Analysis 

3.1. 1.3 Millimeter Emission 

Figure [T] shows a contour image of the 1.3 millimeter emission overlaid on a Hubble 
Space Telescope/ AC S coronographic image of optical scattered light (F606W filter) from 



Graham et al.l (120071 ). The synthesized beam size for this 1.3 millimeter image is 3'.'3 x 2'.'6 



(33 X 26 AU), position angle —12°, obtained with natural weighting and an elliptical Gaussian 
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taper I'.'S x I'/O (FWHM) oriented east-west to make the beam shape more circular. The 
rms noise in this image is 0.40 mJy beam~^, and the individual peaks have a signal-to- 
noise > 7. The star symbol marks the stellar position corrected for proper motion, offset 
by (3'.'20, — 4'.'25) from the phase center. The image shows a resolved band of 1.3 millimeter 
emission that extends approximately symmetrically from the stellar position to the southeast 
and northwest. By fitting a line to the positions of the two peaks in the image, we estimate 
the position angle of the emission structure to be 130 ± 2 °, in exce llent agreement with the 
disk orientation inferred from scattered light observations (JLiull2004l ) . From its position angle 
and double peaked morphology, we identify the emission structure as a limb-brightened dust 
belt. Note that thermal emission from the stellar photosphere contributes only 0.06 mJy 
at this wavelength, and even the s trongest synchrotr on radio flares from stellar activity 
should remain well below the noise (JKundu et al.l 119871 ). The total flux density obtained by 
integrating over a box that surrounds the emission is 8.5 ± 2 mJy, where the uncertainty is 
estimated conservatively by evaluating boxes of the same size in emission-free regions of the 
image. No other significant features are detected in the field of view. 

An extrapolation of the measurements at shorter wavelengths indicates that missing 
flux due to the spatial filtering by the SMA is not a sigri ificant issue for these observatio ns. 
The spectral index between 350 yum (jChen et al.l |2005| ) and 850 /xm (JLiu et al.l |2004| ) is 



l-7lo'5, which predicts 7.1^^2 mJy at 1.3 millimeters, in good agreement with the SMA 
measurement. This shallow index is consistent , within the uncerta inties, with values near 
2.0 determined for a sample of debris disks by iGaspar et al.l (1201 ll ). notably for stars later 
than F-type where dust temperatures are so low that the Rayleigh- Jeans approximation in 
this part of the spectrum is invalid. We conclude that the compact structure detected by 
the SMA accounts for all of the 1.3 millimeter emission in the AU Mic system. 



3.2. A Simple Disk Model 



To characterize the millimeter e mission structu r e, we use a simple parametric disk model, 
following closely the method used by IWilner et al.l (120111 ) for analyzing similar observations 
of the /3 Pic disk. We assume the emission arises from a geometrically thin, axisymmetric belt 
with a radial surface brightness profile given by / oc f{r)r~'^. This functional form is intended 
to capture the essence of low optical depth emission from a disk with a spatially invariant 
dust emissivity, surface density profile /(r), and temperature profile that falls off with radius 
as a power-law. We fix the power-law index q = 0.5 to approximate radiative heating from 
the central star, and we consider two shapes for /(r): (1) an annulus with power law slope, 
/(r) oc r'P for R ± AR/2, p E 0,1, and (2) a Gaussian, /(r) oc exp [-((r - R)/V2ARf]. 
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In each case, the model is described by three parameters: a center radius R, width AR, 
and total flux density F = J IdQ. We fix the dis k inclination at 8 9?5 and position angle 



at 130°, as determined from scattered light data (JKrist et al.ll2005l ). and place the origin 



within O'.'l of the nominal stellar position (within the SMA astrometric accuracy). Note that 



the v ertical thickness of the millimeter grain population is expected to be < 0'.'25 (JThebault 



20091 ). negligible compared to the resolution of the data. The utility of this simple modeling 
approach is to provide basic constraints while avoiding the multitudinous assumptions about 
grain properties, dynamics, and radiative transfer required in more sophisticated treatments. 

To estimate the model parameters and their uncertainties, we calculate three-dimensional 
grids spanning appropriate ranges of the parameter values. For the Gaussian model, the grid 
covers 20 < i? < 47 AU and 1 < AR < 40 AU in 1 AU steps, and 5 < F < 11 mJy in 
0.2 mJy steps. For each position in a model grid, we compute a set of synthetic SMA vis- 
ibilities and compare directly to the data with a x^ value (the sum of real and imaginary 
components over all spatial frequencies) that uses natural weights modified by the Gaussian 
taper used for the image in Figure [TJ The rightmost panel of Figure [2] shows the x^ sur- 
face in the {R, AR) plane for the Gaussian model after marginalizing over F. The best-fit 
model is marked by a cross, and the contours delineate (marginalized) la intervals in Ax^- 
The best-fit center radius is i? = 36l]^g AU, with corresponding width AR = lOlg^ AU (or 
FWHM = 2V21n2Ai? = 23+^^ AU) and fiux density F = 8.0 ± 1.2 mJy. The three lefthand 
panels of Figure |2] show the original 1.3 millimeter image from Figure [11 the Gaussian model 
image obtained using the best-fit parameter values and the same visibility weighting scheme, 
and the imaged residuals (where data— model subtraction is conducted on the visibilities). 
This comparison shows clearly that the best-fit model reproduces the main features of the 
image and leaves no systematic residuals. The power-law models for /(r) give similar best-fit 
parameter values: for p = 0, we find i? = 35 ± 6 AU and A_R = SStH AU (where the upper 
limit is set by the computational grid boundary). For p = I, the best-fit locus is shifted to 
roughly 10% larger values of R, a manifestation of the degeneracy between the radial gradi- 
ent and emission extent. All of these best-fit models reproduce the observed SMA visibilities 
equally well (i.e., have effectively the same minimum x^ values). 

The modest signal-to-noise of the observations preclude placing tighter constraints on 
the model parameters. Nonetheless, it is reassuring that consistent results are obtained from 
the various assumptions about the emission morphology. Each gives a similar center radius, 
as well as an effectively c leared central cav ity with size compatible with that inferred from 



scattered light (~ 12 AU iKrist et al.l 120051 ) and the lack of mid-infrared emission (~ 17 AU 
Liu et al.ll2004l ). The best-fit Gaussian model reproduces better the contrast of the maxima 
to the extended emission of the disk in the image plane, but there is no clear preference for 
any of the functional forms. Of course, none of these ad hoc models are perfect, and more 
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sensitive data will be needed to constrain further the details of the emission distribution. 



Discussion 



We have spatially resolved the thermal emission at 1.3 millimeters from the AU Mic 
debris disk, thereby revealing the distribution of its millimeter grain population. These large 
grains are expected to have dynamics similar to the unseen dust-producing planetesimals, 
and the emission at this long wavelength provides a direct link to their spatial distribution. 

The overall structure of the AU Mic millimeter emission is reproduced well by a belt sur- 
rounding the star centered near a radius of 35 AU. The inferred location of this emission belt 
coincides with the region over which the scattered light brightness profile of the disk steep- 
ens markedly, the feature that provided i mpetus for developi ng the "birth ring" model. Like 
AU Mic's more massive sister star /3 Pic ( IWilner et al.ll201ll ). the resolved multi-wavelength 
data are consistent with a scenario where destructive collisions of planetesimals within the 
belt create grains with a spectrum of sizes, and the effects of size-dependent dust dynamics 
generate the millimeter emission and spatially extended scattered light nebula. For P Pic, 
the intense stellar radiation pressure can account for the large halo of small grains in the 
outer disk. But for AU Mic, which is about two orders of magnitude less luminous, the much 
weaker stellar radiation field is apparently augmen ted by a stellar wind to expel small grains 



weaker stellar radiation iieid is apparently augmen ts 
fetrubbe fc ChianeJ bood : lAueereau fc Beustlbooeh . 



Though the fractional radial extent of the AU Mic millimeter emission is not strongly 
constrained, the best-fit models suggest that it could be broad, AR/R ~ 1. For f3 Pic, the 
millimeter emission is characterized by AR/R < 0.5, with the upper end favored for the width 
of the underlying dust -producing planetesimal belt from detailed models of the knee of the 
scattered light profile ( lHahnll2010l ). Such broad belts are not unusu al features of debris disks, 
in particular among the few resolved at millimeter wavelengths (JHughes et al.ll201l[ ). An 
inner cavity devoid of dust is plausibly maintained by the combination of grinding collisions 
in the belt and clearing by stellar radiation forces, though a tenable alternative is that the 
inner regions are swept clear by a planet. Planets seem to be the most viable explanation for 
the millimeter cavities common to t he gas-rich transition disks around younger, pre-main- 
sequence stars (j Andrews et al.l 120111). While a plane t has been imaged to orbit within the 
cavity in the disk around /3 Pic ( jLagrange et al.ll2010l ). direct evidence for p lanets in the AU 



Mic system remains elusiv e, both from high contrast imaging observa tions ( iMasciadri et al. 



20051 : iMetchev et al.ll2005l ) and from transit searches ( iHebb et al.ll2007l ). The outer boundary 
of the emission belts may refiect the outer extent of successful planetesimal formation, or 
it may simply mark the edge of the region where the planetesimal stirring mechanism has 
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been s uccessful at initiating dust-producing collisions, as in the models of iKenyon &: Bromley 



too4 \ 



These new millimeter observations of AU Mic show the global structure of large grains 
in its debris disk, but they do not have sufficient resolution and sensitivity to reveal any 
significant departures from axisymmetry or substructure that might result from gravita- 
tional interactions with unseen planets. The higher resolution images of scattered light 
clearly show several radial and vertica l inhomogeneities at su barcsecond scales (features A- 



E, see iLiul 12004 : iMetchev et al.l 120051 : Fitzgerald et al.l 120071 ). If these substructures have 



counterparts at millimeter wavelengths, then future observations with the Atacama Large 
Millimeter/Submillimeter Array (ALMA) may be able to detect them, pinpoint the location 
of parent bodies, and constrain their origins. 

We thank James Graham for providing the Hubble Space Telescope image in Figure 1. A. 
M. H. is supported by a fellowship from the Miller Institute for Basic Research in Science. 
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Fig. 1. — SMA image of the 1.3 millimeter continuum emission fr om AU Mic, overlaid on an 
image of optical scattered light from the Hubble Space Telescope (JGraham et al.l 120071 ). The 
contour levels are —2,2,4,6 x 0.40 mJy (the rms noise level); negative contours are dotted. 
The ellipse in the lower left corner represents the 3'.'3 x 2'.'6 (FWHM) synthesized beam size. 
The star symbol indicates the location of the stellar photosphere. 



This preprint was prepared with the AAS E^TJtjX macros v5.2. 
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Fig. 2. — Lejt panels: SMA image of the 1.3 millimeter emission from AU Mic together with 
images of the best-fit axisymmetric belt model and residuals (see text). The contour levels 
and beam size are the same as in Figure [H The dashed line indicates the 130° position angle 
of the scattered light disk. Right panel: The x^ surface for the model emission belt center 
and width parameters, with contours at 1, 2, 3a. The cross marks the best-fit model. 



Table 1. Submillimeter Array Observations of AU Mic 



Observation 


2011 Jul 31 


2011 Sep 26 


2011 Sep 27 


2011 Oct 25 


Array Configuration 
Number of Antennas 
Baseline Lengths (m) 
LO frequency (GHz) 
225 GHz atm. opacity*" 


Extended 

8 

10-189 

235.6 

0.10 


Compact 

6 

8-68 

225.4 

0.10-0.06 


Compact 

6 

8-68 

225.4 

0.06 


Compact 

7 

8-68 

235.6 

0.08 



"Measured at the nearby Caltech Submillimeter Observatory. 



